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We report a new approach to the thermal conductivity manipulation - substrate coupling. 
Generally, the phonon scattering with substrates can decrease the thermal conductivity, as observed 
in recent experiments. However, we find that at certain regions, the coupling to substrates can 
increase the thermal conductivity due to a reduction of anharmonic phonon scattering induced 
by shift of the phonon band to the low wave vector. In this way, the thermal conductivity can 
be efficiently manipulated via coupling to different substrates, without changing or destroying the 
material structures. This idea is demonstrated by calculating the thermal conductivity of modi- 
fied double- walled carbon nanotubes and also by the ice nanotubes coupled within carbon nanotubes. 
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With the shrinkage of electronic devices to the 
nanoscale[ll-[3] and the revival of thermoelectrics0-0], 
thermal transport property of nanomaterials has become 
more and more significant. So far, great efforts have 
been done on finding/synthesizing new nanomaterials 
with particularly high/low thermal conductivities. 0413 
Recently, some efforts have taken headway on the sub- 
ject of manipulating the thermal conductivity via doping, 
adsorbing, or generating defects. 
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These processes 

can only reduce but hardly enhance the thermal con- 
ductivity, since the phonon scattering in a crystal lattice 
usually be less than that with the doping or defects. Fur- 
thermore, most of these treatments largely destroy the 
structure and thus the corresponding properties of nano- 
materials, which make them inapplicable in the nanodc- 
vices. 

Very recently, the effects of coupling with substrates on 
thermal conductivity were reported (16H18J. In general, 
when a conductive material is coupled with substrate, its 
thermal conductivity is expected to be decreased owning 
to the additional phonon scattering with the substrate as 
observed in the recent experiments. [l6l. [l7| In this letter, 
we find that, at certain regions, the effect of phonon scat- 
tering can be suppressed and the thermal conductivity of 
nanomaterials can be surprisingly increased due to the 
coupling induced shift of phonon band to the low wave 
vector. Based on this finding, we propose a new approach 
to thermal conductivity manipulation-coupling with dif- 
ferent substrates. This approach naturally extends the 
capability of conventional treatments on thermal con- 
ductivity without destroying the structures of materials. 
Since in the production of nanodevices, the conductive 
nanomaterials are always placed on certain substrates, 
our approach has great potential for advancing the per- 
formance of nanoelectronic and thermoelectric devices. 
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FIG. 1: Schematic configuration of our approach to manipu- 
late thermal conductivity. Chain- 1 and chain-2 represent the 
conductive material and the substrate, respectively. The two 
particles on the left and right ends of chain- 1 are put into 
contact with thermostats of temperature Th and Tl, respec- 
tively. 



In order to demonstrate this approach, we start with 
a coupled Fermi-Pasta-Ulam(FPU) chainJH, H(| model. 
Then we use two examples, thermal conductivity of both 
modified double- walled carbon nanotubes (DWNTs) and 
coupled ice nanotubes (Ice-NTs) [2l[ , to demonstrate the 
applicability in real systems. 

Fig. 1 shows the schematic configuration of our ap- 
proach, which is illustrated by a coupled atom chain 
model. The upside chain (chain-1) represents the con- 
ductive material and the underside chain (chain-2) rep- 
resents the substrate. The two ends of chain-1 are con- 
tacted with the thermostats with temperature Th and 
Tl, respectively, while, chain-2 is free of thermostat con- 
tact. Each node of chain-1, except for the two ends con- 
tacted with thermostat, is coupled to the corresponding 
underside node of chain-2. This ladder-like construction 
corresponds well to the real systems, where the conduc- 
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FIG. 2: Heat flux of chain-1 and chain-2 with variation of 
resonance angle $ in the coupling system for three cases: (a) 
=1, ni2 is varied; (b) m2=l, &2 is varied; (c) both 7712 
and k2 are varied but keeping 7712 =fe- The chain length iV 
is 50, and Th—0.3, Tl=0.2. The solid squares, triangles, and 
circles represent the heat flux of chain-1 for cases (a), (b), 
and (c), respectively; the corresponding open ones represent 
heat flux of chain-2 for the three cases. The short dashed 
line shows the heat flux of isolated chain-1. Chain- l's heat 
flux monotonously decreases with increasing, which can be 
either larger or smaller than that of the isolated one depending 
on the specific * value. The inset shows (m2, fe) dependence 
of the resonance angle for the three cases. The resonance 
angle is very sensitive to the (m.2, fe) parameters. 



tive materials are always fixed on some substrates. 

Two coupled FPU chains are firstly considered to rep- 
resent the approach in details, with all the Hamiltonian 
parameters being in the reduced units. [19l. [20( In the 
model, the atomic mass and spring constant of chain- 
1 (mi, fci) are fixed to be 1.0, the anharmonic coefficient 
of chain-1 and chain-2 (/3i, (3 2 ) and the coupling strength 
k c are all set to be 0.5. The atomic mass and spring 
constant of chain-2 (m 2 ,k 2 ) are variables to manipulate 
thermal conductivity of chain-1 (see Supporting Informa- 
tion, SI.I). 

The nonequilibrium molecular dynamics fNEMD) [l9l. 
I20I [22j method is used to calculate the heat flux of chain- 
1 and chain-2 in the coupling system. It is found that 
the phonon resonance between chain-1 and chain-2 plays 
an important role on the thermal conductivity. Here 
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to describe 



we use the resonance angle 

the phonon resonance strength, where Ai and A2 are the 
phonon amplitudes of chain-1 and chain-2 after coupling 
(see Supporting Information, SI. II) . It is obvious that the 
resonance strength becomes a maximum (minimum) at 
*=f(* = 0). 

Fig. 2 shows heat flux J of chain-1 and chain-2 with 
variation of the resonance angle. Since VP is not uniquely 
determined by m 2 or k 2 , many (iri2, k 2 ) sets can be cor- 
responding to a $ value, we mainly consider three cases: 
(a) k 2 =1, m 2 is varied; (b) m 2 —l, k 2 is varied; (c) both 



7712 and k 2 are varied but keeping m 2 =k 2 . For all the 
three cases, the heat flux of chain-1 monotonously de- 
creases with \I> increasing, while the heat flux of chain-2 
has a contrary behavior. Particularly, when the phonon 
resonance becomes strong enough (with 1 F nearby j), 
heat flux of chain-1 gets much smaller than that of the 
isolated one (J = 0.034), suggesting a substantial reduc- 
tion of thermal conductivity [23J; when the phonon res- 
onance becomes small enough, however, the heat flux 
of chain-1 gets obviously larger than that of the iso- 
lated one, showing an increment of thermal conductiv- 
ity, which has not been explored so far. This interesting 
phenomenon shows that the thermal conductivity of con- 
ductive material can be efficiently manipulated through 
the substrate. 

We note that the resonance angle is very sensitive to 
the (m 2 , k 2 ) parameters (Inset of Fig. 2). This suggests 
that one can easily find a proper substrate to manipulate 
the thermal conductivity in the application. For example, 
to decrease the thermal conductivity, the substrate needs 
to have similar atomic mass and spring constant with the 
conductive material; to increase the thermal conductiv- 
ity, the atomic mass and/or spring constant of substrate 
only needs to be several times larger than that of conduc- 
tive material. Moreover, the heat-flux curves of cases (a), 
(b), and (c) are very close to each other, indicating that 
the thermal conductivity of chain-1 is insensitive to the 
specific type of {m 2 , k 2 ) compositions. Thus there would 
be rich choices of substrate candidates for the thermal 
conductivity manipulation. Consequently, the present 
approach can be easily realized in the nanotechnologies, 
which provides a clear direction for designing advanced 
nanoelectronic and thermoelectric devices. 

The reduction of thermal conductivity after coupling 
can be understood from the phonon-resonance effect, 
which induces strong phonon scattering and thus reduces 
the thermal conductivity. [HI l24T - |27j To understand the 
increment of thermal conductivity, however, we need to 
invoke phonon band theory. Here we consider an extreme 
condition, where the resonance angle is zero and the in- 
crement of thermal conductivity is maximum. From the 
coupling-harmonic-oscillator (CHO) model, the phonon 
dispersion of chain-1 in the coupling system can be writ- 
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ten as uj(q) 



' 4 sin 2 (|g) + k c . Compared with 



that of the isolated case uj(q) = 2y ^ ■ sin(ig), the 

phonon band of chain-1 has an obvious upshift after cou- 
pling, the magnitude of which is proportional to \fk~ c (In- 
set of Fig. 3). 

Based on the phonon dispersion, the phonon energy 
density (PED) distributions of chain-1 with wave vector 
q can be further calculated (see Supporting Information, 
SI. III). Compared with that of isolated case, the PED 
peak has an obvious shift to the small- q direction own- 
ing to the upshift of phonon band after coupling (Fig. 
3), suggesting more energy has been carried and trans- 
ported by the small-g phonons. The phonons' scattering 
power is proportional to the phase difference between 
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FIG. 3: PED distributions with the wave vector q of chain-1 
in both the isolated and the coupled (with \I/ = 0) cases at 
T=0.25. The PED peak has an obvious shift to the small-g 
direction after coupling. The inset shows the corresponding 
phonon dispersions. 
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FIG. 4: (a) Temperature dependence of heat flux of isolated 
chain-1 (iV=50). The phonon scattering effect has become 
dominant at our simulation temperature (T=0.25). The inset 
shows chain length N dependence of total heat flux JN of 
chain-1 in both the isolated and the coupled (with ^ = 0) 
cases. The larger fitting exponent implies that the increment 
of thermal conductivity by the coupling gets more and more 
distinguished with the chain length increasing. 



different atoms, which is hence proportional to the q 
value. [Ill More small-g phonons being responsible for the 
heat transport corresponds to smaller scattering power, 
thus larger phonon mean free paths (PMFPs) in chain- 
1, which increases the thermal conductivity (positive ef- 
fect). On the other hand, the upshift of the phonon band 
also results in less phonons being excited for the heat 
transport, which in turn reduces the thermal conductiv- 
ity (negative effect). The thermal conductivity varia- 
tion induced by the phonon band upshift is attributed 
to such two effects competing with each other (phonon- 
band- upshift effect). 



We have also calculated the temperature dependence 
of heat flux of isolated chain-1 (Fig. 4) and found the 
heat flux monotonically decreases with the temperature 
increasing from T=0.15 to 0.65. It is known that in- 
creasing temperature has two competitive effects on the 
thermal conductivity: (i) It excites more high-frequency 
phonons that enhance the thermal conductivity, (ii) 
It increases phonon-phonon scattering that reduces the 
thermal conductivity. [HI, [2(| Fig. 4 shows that effect 
(ii) has become dominant at our simulation temperature 
(T=0.25), indicating that the upshift of phonon band, 
which has similar effect with that of decreasing temper- 
ature, has more positive effect than the negative effect 
on the thermal conductivity. Hence the thermal conduc- 
tivity can be increased due to the phonon-band-upshift 
effect. Moreover, inset of Fig. 4 shows the chain length 
dependence of total heat flux(JiV) of chain-1 in both the 
coupled (with >F = 0) and isolated cases. As is shown, 
JN of the isolated chain-1 diverges as JN oc TV ' 4 , consis- 
tent with the results in previous reports. [li| Ho| While, 
for the coupled case, JN oc N 0A8 . The larger fitting 
exponent implies that the increment of thermal conduc- 
tivity by the coupling gets more and more distinguished 
with the chain length increasing. 

Consequently, in the coupling system, the thermal con- 
ductivity variation is owed to both the phonon-band- 
upshift effect and the phonon-resonance effect that com- 
pete with each other. From the results above (Fig. 2), 
the phonon-resonance effect can be easily manipulated by 
changing the atomic mass and/or spring constant of the 
substrate. Thus we can efficiently manipulate the ther- 
mal conductivity of conductive material through coupling 
it to different substrates. 

As a demonstration of this approach in the real sys- 
tems, we have calculated the thermal conductivity of a 
(7,0) carbon nanotube (CNT) coupled within (14,0) CNT 
as a substrate((7,0)@(14,0) DWNT, see Supporting In- 
formation, SI.IV). We set the atomic mass of (14,0) CNT 
(Month) as parameter to see how the thermal conduc- 
tivity of (7,0) CNT changes with the resonance angle. 
In the coupling system, the resonance angle 'F and thus 
the thermal conductivity can be efficiently manipulated 
through changing Month- As shown in Fig. 5(a), ther- 
mal conductivity of (7,0) CNT obviously increases with 
Month increasing from one to several times that of 
atomic mass of (7,0) CNT (carbon-atom mass, Mcntt) 
owing to the reduction of phonon-resonance effect. Com- 
pared with that of the isolated case, the thermal conduc- 
tivity of (7,0) CNT can be either substantially decreased 
or increased depending on the value of Month, consis- 
tent with the results of coupled FPU chain model, which 
can be well understood form the coupling mechanism dis- 
cussed above. When Month = Mcnt7 ~ f ), the 
phonon-resonance effect that decreases the thermal con- 
ductivity is dominant and the thermal conductivity is 
decreased by coupling. With Month increasing, the 
phonon-resonance effect becomes more minor. When 
Month > 2M C nt7 (* < ^)[29|. the phonon-band- 
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FIG. 5: (a) Thermal conductivity of (7,0) CNT coupled 
within (14,0) CNT, with variation of (14,0) CNT's atomic 
mass Mqntu [Monti is kept to be carbon-atom mass). The 
short dashed line shows the thermal conductivity value of 
isolated (7,0) CNT. The thermal conductivity of (7,0) CNT 
can be either substantially decreased or increased depending 
on the value of Month- The inset shows the structure of 
(7,0)@(14,0) DWNT. (b) Thermal conductivity of Ice-NTs. 
Thermal conductivity of Ice-NTs has an obvious increment 
after coupling. The inset shows optimized structures of ice- 
NTs coupled within CNTs. 



ity of Ice- NTs is expected after coupling with CNTs. Fig. 
5(b) shows the calculated thermal conductivity of various 
Ice-NTs both with and without CNT coupling (see Sup- 
porting Information, SI.V). As one can see, the thermal 
conductivity of Ice-NTs has an obvious increment after 
coupling, which is independent on the specific tube types. 
These realistic illustrations further confirm the feasibility 
of our approach for applications. 

In summary, we have proposed a new approach to ma- 
nipulate the thermal conductivity. By coupling with dif- 
ferent substrates, thermal conductivity of the conduc- 
tive nanomatcrial can be either remarkably decreased or 
increased, which can be realized in the device applica- 
tions. To decrease the thermal conductivity, the sub- 
strate needs to have similar atomic mass and spring con- 
stant with the conductive material; to increase the ther- 
mal conductivity, the atomic mass and/or spring con- 
stant of the substrate needs to be several times larger 
than that of the conductive material. Through the illus- 
trations of double-walled carbon nanotubes and coupled 
ice nanotubes, we have further shown that this approach 
is applicable in the real systems. Compared with the 
conventional treatments which only reduce the thermal 
conductivity, our approach can truly realize the thermal 
conductivity manipulation in solid nanomaterials with- 
out destroying their structures. 



upshift effect that increases the thermal conductivity be- 
comes dominant, and thus the thermal conductivity is 
increased by coupling. 

Ice-NT coupled within a CNT can be considered as 
another realistic illustration for our approach, where the 
Ice-NT and CNT correspond to the conductive mate- 
rial and the substrate. Since the CNTs have much 
larger spring constant than the Ice-NTs (from Supporting 
Information,^ ~ 0), an increment of thermal conductiv- 
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I. SUPPORTING INFORMATION 

" Manipulating thermal conductivity through substrate 
coupling" , Zhixin Guo, Dier Zhang, and Xin-Gao Gong 
SI. I Hamiltonian of the coupled FPU chains 
The Hamiltonian of the coupled FPU chain model can 
be written as: 



H = H 1 +H 2 +H C 



(1) 



amplitude ratio of chain-2 and chain-1 after coupling can 
be expressed as: 



■ tan(a) (m,2 > mi), 



(4) 



where Ai and A2 are the amplitudes of chain- 
1 and chain-2, respectively. Here tan(a) = 

fee 



^3 c - € [-1,1], andw, 



fcl+fcc f. _ k 2 + k c 
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gle is defined as ^ = 



V m l m 2 : 

respectively. Then the resonance an- 
A2 ^ When * = 0, the 



arctan(^ 

phonon dispersion of chain-1 in the coupling system can 
be written as: 



u(q) 




4sin 2 (-g) + k c 



(5) 



Thus the phonon band of chain-1 has an obvious upshift 
after coupling, the magnitude of which is proportional to 

SI. Ill Calculation of PED distribution with wave vec- 
tor 
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PED distributions with phonon frequency u) of iso- 
cfiain-1 at T=0.25. 



where H\,H2, and H c are the Hamiltonian of chain-1, 
chain-2, and the coupling term, respectively. x l n , p l n 
represent the displacement from the equilibrium position, 
the momentum of the ith particle; and m n , k n , Pn are 
the mass, spring constant and the anharmonic coefficient 
of chain-n (n = 1,2). 

SI. II Definition of resonance angle ^ 

Using the coupling-harmonic-oscillator (CHO) model, 
we can calculate the phonon resonance strength. The 



The q dependence of PED has an expression of 
u(q) — hcu(q) x r)[u!(q), T],with r)[u)(q),T] representing 
the phonon distribution function at temperature T for 
a certain frequency uj. Since the phonon distribution in 
a heat-transport system is mainly determined by the heat 
source that supplies phonons. We suppose the conduc- 
tive chain is coupled to an ideal heat source that supplies 
phonons like the black-body, which can be represented by 
the Nose-Hoover thermostat P, 0]. Then rj(uj,T) can be 
written as: 



87TW 2 



1 



„3 e hcu/k b T _ 1 ' 



(6) 
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Using equation (6), we also calculated the PED distri- 
butions with phonon frequency uj of the isolated chain- 1. 
As is shown in Fig. 1, the result is consistent with that 
calculated by Li et al. using the numerical method, which 
confirms the reliability of our consideration. [3| 

SI. IV NEMD simulation for the double- walled carbon 
nanotube 

In the simulation, the Tersoff potential was used to 
the describe the carbon-carbon inter action and the 
coupling interaction between two CNTs was described 
by the LJ potential, [5] which was truncated at 10 A by a 
switching function. The simulated (7,0)@(14,0) DWNT 
has a length of 121 A, containing 2352 atoms. The wall- 
thickness of CNT was chosen to be 1.44 A for the calcula- 
tion of cross-section area. 0, 0] Fixed boundary condition 
was applied along the axial direction of CNTs, where the 
outmost two layers of each head were fixed. @, H| Then 
two layers of each end of (7,0) CNT were put into contact 
with the Nose-Hoover thermostat with temperatures 310 
and 290 K, respectively, while the (14,0) CNT was free 
of thermostat contact. To integrate equations of motion, 
the velocity Verlet method was employed with a fixed 
time step of 1 fs. All results were obtained by averaging 
about 10 ns after a sufficient long transient time (10 ns) 



to set up a nonequilibrium stationary state. 

SI.V NEMD simulation for the Ice-NTs 

Water-water intermolecular interaction was described 
by the TIP4P[9] potential and carbon-carbon interaction 
was described by the Tersoff[4| potential. As for the cou- 
pling term, the CNT- water interaction was described by a 
carbon-oxygen L J potential. fic| All the pair interactions 
were truncated at 10 A by a switching function. The sim- 
ulation box has a length of 121 A and the total number 
of water molecules inside is 44 xn for the (n,0) Ice- NT, 
where n = 4, 5, and 6. The Ice-NTs' wall-thickness was 
chosen to be 2.75 A for the calculation of cross-section 
area[ll|. Fixed boundary condition was applied along 
the axial direction of Ice-NT/CNT, where the outmost 
two layers of each head were fixed. @, Q Then two lay- 
ers of each end of Ice-NT were put into contact with the 
Nose-Hoover thermostat with temperatures 110 and 90 
K, respectively, while the CNT was free of thermostat 
contact. To integrate equations of motion, the velocity 
Verlet method was employed with a fixed time step of 1 
fs. All results were obtained by averaging about 5 ns 
after a sufficient long transient time (5 ns) to set up a 
nonequilibrium stationary state. 



[1] Nose, S. J. Chem. Phys. 1984, 81, 511-519. 
[2] Hoover, W. G. Phys. Rev. A 1985, 31, 1695-1697. 
[3] Li, B., Lan, J.; Wang, L. Phys. Rev. Lett. 2005, 95, 
104302. 

[4] Lersoff, J. Phys. Rev. B 1989, 39, 5566-5568. 

[5] Girifalco, L. A.; Hodak, M.; Lee, R. S. Phys. Rev. B 

2000, 62, 13104-13110. 
[6] Zhang, G.; Li, B. J. Chem. Phys. 2005, 123, 114714. 
[7] Guo, Z.; Zhang, D.; Gong, X. G. Appl. Phys. Lett. 2009, 

95, 163103. 



[8] Wu, G.; Li, B. Phys. Rev. B 2007, 76, 085424. 

[9] Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; 

Impey, R. W.; Klein, M. L. J. Chem. Phys. 1983, 79, 

926-935. 

[10] Meng, L.; Li, Q.; Shuai, Z. J. Chem. Phys. 2008, 128, 
134703. 

[11] Guo, Z.; Zhang, D.; Zhai, Y.; Gong, X. G. Nanotechnol- 
ogy 2010, 21, 285706. 



